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ABSTRACT 

This dynamical model for the 28 galaxies with distances less than 1.5 Mpc, 
and not apparently tight satellites, is constrained by the initial condition that 
peculiar velocities at high redshift are small and growing in accordance with the 
standard cosmology. The solution is a satisfactory fit to most of the measured 
redshifts, distances, and proper motions, with some interesting exceptions that 
call for further investigation. The model predicts Milky Way rotation speed 
256 km s~^, consistent with Reid et al. (2009a). Ten Local Group galaxies 
emanate from low supergalactic latitude and supergalactic longitude ~ 70°, per- 
haps as remnants from failed assembly of a larger galaxy. NGC 6822 passes close 
to the Milky Way at redshift z ~ 0.27, in an orbit similar to the Magellanic 
Clouds. Leo 1 has heliocentric angular velocity 0.33 mas yr^^, perhaps measur- 
able by the mean stellar motion, and 15 galaxies have proper motions greater 
than 0.05 mas yr~^, measurable for any with masers. 

1. Introduction 

Advances in the measurements of distances and proper motions of nearby galaxies arc 
producing a rich data set and the opportunity to explore dynamics in the local universe in 
much greater detail than was possible a decade ago (as in Peebles, Phelps, Shaya, & Tully 
2001). This is a report of progress in the analysis of what happened within the Local Group 
(LG) based on the redshift and distance data in the Local Universe Catalog, also a work 
in progress (RBT and EJS), and on the measured proper motions of three Local Group 
(LG) members. The dynamical model includes orbits for the Milky Way (MW) and the 27 
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galaxies with measured distances less than 1.5 Mpc, after excluding the Small Magellanic 
Cloud, whose orbit may be entangled with the Large Magellanic Cloud (LMC), and excluding 
the galaxies that are judged to be too close to MW or M31 for a feasible analysis of where 
they were at high redshift. The circular velocity Vc of MW at the Solar circle also is a 
parameter, mainly for the conversion between heliocentric and galacto centric velocities (as 
noted by Shattow & Loeb 2009). 

The starting condition from cosmology is that the peculiar velocities of the matter at 
high redshift are generally smaller than now and increasing in accordance with the linear 
perturbation theory of growing departures from a homogeneous mass distribution of an ideal 
pressureless fluid. The resulting numerical solutions for the orbits of LMC and M33 at 
low redshift are similar to solutions that ignore this initial condition (Besla et al. 2007; 
Kallivayalil et al. 2009; Putman et al. 2009), because the present positions and velocities 
are similar. Orbits computed with and without the cosmological initial condition diverge at 
high redshift, however, in analogy to the difference at high redshift between a growing mode 
and a mixture of growing and decaying modes in perturbation theory. LG analogs derived 
in numerical simulations (as in Boylan-Kolchin, Besla, & Hernquist 2010; Libeskind ct al. 
2011) also satisfy the cosmological initial condition, and are an important complement to 
this explicit analysis of what is observed in the Local Group. 



2. The Model 

2.1. Approximations 

This analysis depends on three major approximations. The flrst is the representation 
of each LG galaxy by a point particle back to the starting time of the computation, at 
expansion factor 1 -\- zi — 10. The physical picture is that each particle is meant to track 
the mean motion of the fragments that are gathering and merging into a protogalaxy, and 
to represent the gravitational attraction of all these fragments on neighboring protogalaxies. 

The second major approximation is the initial condition at Zi. In linear perturbation 
theory at high redshift, z, the growing density contrast in an ideal pressureless fluid scales as 
5p/p oa 1/(1 + z), and the peculiar velocity v and acceleration g are related to the expansion 
time t from 2; — >■ oo by 

v^gt. (1) 

We apply this initial condition to the particle model. An equivalent condition is that the 
comoving coordinate displacements from the particle positions ait — scale as x{t) — a;(0) oc 
a{t). This is an approximation, among other reasons because it does not take account of the 
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nonlinear development of structure on the scale of galaxies at Zi. The point is that it fits 
the idea that peculiar velocities of protogalaxies at high redshift are small and growing as 
structure grows. 

The third approximation is the use of four external actors for a phenomenological de- 
scription of the effect of matter external to LG on the motions of the LG galaxies. The 
actors represent in first approximation the neighboring mass concentrations in the Sculptor 
Group and the Maffei and Centaurus systems, but the masses, positions and redshifts of the 
actors are allowed to float to improve the model fit to the measurements of the LG galaxies. 
The approximation would fail if a large external mass that now has little effect on LG was 
a serious influence in the past, when these masses were much closer to LG. This need not 
be a problem, of course, because LG was more compact too. And we have a meaningful in- 
dication that this and the other two approximations yield a useful approach to reality, from 
the generally reasonable fit of the model to considerably more measurements than there are 
relevant adjustable parameters. But an important goal for a more complete analysis is to 
check the third approximation by taking explicit account of the dynamics and gravitational 
effects of the observed external mass concentrations. 



2.2. Numerical Method 

In the form applied here the numerical action method (NAM) of dealing with the mixed 
boundary conditions of given present positions and the initial condition in equation ([T]) 
represents particle orbits by positions at discrete time steps (Peebles 1989, 1995; Peebles, 
Phelps, Shaya, &; TuUy 2001). It produces a solution to the equation of motion in leapfrog 
approximation that is reached by iterated joint shifts of all particle positions from random 
trial orbits toward a stationary point of the action in the direction indicated by the first and 
second derivatives of the action. The method is described in Appendix A. 

The mixed boundary conditions allow multiple solutions, and there is no guarantee that 
our adoption of the best fit to the constraints chooses the most realistic solution. Thus the 
numerical solution presented here certainly may have some wrong orbits. Wrong orbits for the 
smaller galaxies that behave as test particles may be discovered from improved measurements 
of positions and velocities, or perhaps from signatures of disturbances by close passages of 
more massive galaxies. The orbits of more massive galaxies are more tightly constrained by 
their gravitational effects on other orbits, but that is not a guarantee of uniqueness. Also 
to be born in mind is that our measure of fit to the constraints as a function of the 
parameters has many local minima. Here again there is no guarantee that we have arrived 
at the minimum corresponding to the most accurate solution. The straightforward way to 
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explore both ambiguities — multiple solutions and multiple minima of the measure of fit to 
the data — is to repeat construction of the numerical solution starting from many different 
random orbits. This was only very partially done for the solution presented here, in the 
following procedure. 

The construction commenced with dynamical solutions for eight galaxies: MW and 
M31, the three galaxies with measured proper motions — LMC, M33, and ICIO — and 
with NGC 3109, NGC 300, and the Maffei system serving in effect as external actors. The 
small particle number made it practical to find many solutions starting from random or- 
bits (quadratic curves with random coefficients that allow the orbits to arc over several 
megaparsecs) and measured parameters randomly chosen within their assigned standard de- 
viations. In the few solutions with promising fits to the measurements all the free parameters 
were adjusted to minimize a measure of fit. The solution that ended up with the smallest 
was the starting point for the addition of galaxies one at a time. Each addition started 
with a large set of trial random orbits for the new galaxy, and each trial orbit relaxed to a 
solution to the equation of motion for all the galaxies. In the solution with the smallest x^ 
all parameters for all particles were adjusted to minimize x^- In some cases this cycle was 
repeated many times in a search for the best fit to the measurements. 

The value of x^ as a function of the masses, present positions, and Vc has discontinuities 
where a slight change in a parameter causes the solution to jump to quite different orbits, 
generally with much larger x^- To deal with this the sohition was stored before each trial 
adjustment of a parameter, and the solution and parameter restored if the adjustment crossed 
a discontinuity indicated by a distinctively large number of iterations to a solution to the 
equation of motion. A parameter shift that increased x^, or encountered a discontinuity, was 
multiplied by 0.75 and the sign changed for the next round of parameter shifts. A shift that 
lowered x^ was multiplied by 1.25 for the next iteration. 

It would be informative to repeat this whole procedure with different random orbits and 
a different sequence of addition of the galaxies. The practical problem with the numerical 
method used here is that the computation time to add a galaxy to the solution scales as 
the fourth power of the number, n^, of galaxies (varying as for the joint relaxation to a 
solution and as Up to explore how the parameters of all the galaxies in the sohition are to be 
adjusted to minimize the x^ measure of fit). Also, with increasing rip the discontinuities in 
X^ appear more frequently, and discovering them requires more iterations in the relaxation 
to a stationary point of the action. The result was computation time of about a week (on an 
iMac quad) to add the last two galaxies. Sextans A and B, in the solution presented here. 
Use of the present numerical method to explore the many local minima of x^ and the many 
choices of multiple solutions certainly would not be undertaken lightly. But the numerical 
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method is a work in progress, along with the measurements of distances and velocities. 

2.3. Parameters 

Since the model is at best a useful approximation — galaxies are not point particles — 
it cannot be expected to yield a highly precise fit to the galaxy positions and velocities. The 
strategy for dealing with this along with the measurement uncertainties is to assign what 
seem to be moderately optimistic goals for the differences between model and measurements, 
use these goals as effective standard deviations in a measure of fit, and adjust parameters 
to minimize this x^. 

The 28 LG galaxies listed in Table 1 are all the entries in the Local Universe Catalog 
that are closer than 1.5 Mpc and far enough away from MW and M31 that their orbits 
are judged likely to be computable in NAM. The catalog angular positions agree with the 
NASA/IPAC Extragalactic Database (NED), the redshifts and distances occasionally are 
more recent than the NED entries, and many of the catalog distance uncertainties are quite 
different from the NED entries. 

The effective standard deviations in the LG galaxy redshifts are taken to be 5 km s~^ 
for M31, which plays a central dynamical role, and for LMC, whose proper motion is a key 
constraint, with the looser goal of 10 km s^^ difference between model and measurement 
for all the other LG galaxies. This nominal uncertainty is far larger than the precision of 
some redshift measurements, but we judge it to to be optimistic for our numerical model. 
The assigned effective standard deviations in distances are 6 kpc for LMC, 50 kpc for M31, 
and twice the stated uncertainties in Table 1 for the difference between model and measured 
distances of the rest of the LG galaxies. Again, this is arguably optimistic considering the 
limitations of the model and the possibility that the difficult art of distance measurements 
has allowed some errors well above the estimated uncertainties. Offsets between model and 
measured present positions perpendicular to the line of sight may be real — a galaxy of 
stars may not be centered on its dark matter halo — or perhaps more likely a fault of 
the approximate model. The assigned standard deviations in the perpendicular offsets are 

— 2 kpc for LMC, because a larger error causes an objectionable error in angular position 
of this nearby galaxy, — 1 kpc for M31, which partially sets the model orientation, and 
D± — 5 kpc for all the other LG galaxies. The logarithms of the particle masses enter with 
standard deviations equivalent to a factor of ten. This very loose constraint on masses allows 
some indication of how well the model might predict the masses that arc most relevant to the 
dynamics, and which of the galaxies have masses so small, as indicated by the luminosities 
in the Local Universe Catalog, that their masses do not matter. Each orthogonal component 
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Table 1: The Local Group Galaxies 
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of each initial peculiar velocity at 1 + 2;^ = 10 is allowed standard deviation 50 km s~^ in x^. 
This choice is based on the consideration that since the protogalaxies were close to touching 
at 1 + 2;^ = 10 the protogalaxy peculiar motions might be expected to have been comparable 
to motions within protogalaxies. The constraint on initial velocities is needed also to prevent 
NAM from producing an occasional solution with quite unreasonably large Vi. 

The conversion between galactocentric and heliocentric velocities uses the fixed Solar 
velocity components U = 11.1, V = 12.2, W = 7.2 km s^^ relative to the local standard of 
rest (Schonrich, Binncy & Dchncn 2010). The circular velocity of the local standard of rest is 
a parameter to be adjusted. In the computation of we adopt the nominal or catalog value 
Vc — 230 ± 10 km s^^, intermediate between a standard estimate, 220 km s^^, and the larger 
result obtained by Reid et al. (2009a). As discussed in the next section the minimization 
of favors a circular velocity larger than nominal and close to Reid et al. Coordinate 
positions and peculiar velocities are referred to the Priedmann-Lemaitre cosmology with 
flat space sections, Hubble parameter Ho — 70 km s~-^ Mpc~^, matter density parameter 
flm — 0.27, and radiation ignored. 

The value of Vc also figures in the gravitational acceleration produced by MW, which 
is modeled as g — v'^/r at distance r < Vg — Gm/vl, where m is the MW mass, with the 
inverse square law at larger separation. This is relevant for LMC and NGC6822. The mass 
distribution in M31 is modeled the same way, with circular velocity fixed at 250 km s^^, but 
as it happens no model orbits pass within for M31. All other galaxies arc treated as point 
particles. None pass close enough to any galaxy except MW to suggest a serious problem 
with this. 



Table 2: Proper Motions 
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We have constraints on the proper motions from KaUivayahl et al. (2006) and Piatek, 
Pryor & Olszewski (2008) for LMC, Brunthaler et al. (2005) for M33, and Brunthaler et 
al. (2007) for ICIO. Table 2 lists the adopted proper motions as the components ji^ and 
jjLs of the heliocentric angular velocity in the directions of increasing right ascension a and 
increasing declination 5. The adopted angular velocity of LMC is from KaUivayahl et al. 
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(2006). The heliocentric angular velocities of M33 and ICIO are based on measured motions 
of masers that must be corrected for the motions of the masers relative to the host galaxies, 
and the conversion from velocities within the galaxy to angular velocities depends on the 
distance to the galaxy. Distances are adjustable parameters in this analysis, but since the 
effect of changing the distance within a reasonable range is small we adopt the correction at 
distance 809 kpc for M33, with a straight mean of the proper motions from the two masers, 
and distance 760 kpc for the one maser in ICIO. The measurement uncertainties are derived 
from the stated uncertainties in the angular velocities of the masers and in their motions 
within the galaxies. We treat the uncertainties in Table 2 as standard deviations in the 
computation of x^. 



Table 3: External Actors 
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The external actors in Table 3 are nominally the observed mass concentrations in the 
nearby Maffei and Centaurus systems, with the galaxies NGC 55 and NGC 300 assigned to 
represent the still closer Sculptor group. It should be understood that numerical values of 
parameters assigned to these external actors are not to be taken as useful approximations to 
the properties of the real systems: they are meant to provide a phenomenological description 
of the effect of external mass on the LG. For this purpose we allow greater latitudes in 
parameters. Each is assigned effective standard deviation of 50 km s~^ in redshift and a 
factor of 50 in mass. NGC 55 and NGC 300 are assigned effective standard deviations of 
200 kpc in radial distance and 100 kpc in perpendicular offset, and the actors representing 
the Maffei and Centaurus systems are assigned standard deviations of 300 kpc in distance 
and 1 Mpc in perpendicular offset. 



2.4. Numerical Accuracy 

The numerical action solution in this model has positions uniformly spaced in the ex- 
pansion parameter a{t), with rix — 500 steps between oi = 0.1 (meaning the initial redshift 
is Zi — 9) and the present at a„^+i = 1. In the units of the computation (1 Mpc, 100 km s~^. 
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and lO^^mQ) relaxation terminates when the sum of the squares of the action gradients 
(eq. |A2| ) of all particles at all time steps is SOS < 10^^^. The last iteration usually leaves 
the sum at SOS ^ 10~^°. This ensures that the effect on of a small change of a parameter 
signifies a change in the fit to the constraints rather than a change in the numerical error in 
the solution. 



20 
15 
10 

5 



-0.06 -0.03 0.03 0.06 

Fig. 1. — Distribution of fractional differences between the components /c = 1,2,3 of the initial 
LG galaxy peculiar velocities in linear perturbation theory and in the NAM solution. 

The test of initial conditions in Figure [T] compares the Cartesian components of the 
velocities at 1 + = 10 of the 28 LG galaxies in the NAM solution to the linear perturbation 
theory relation v = gt in equation ([T|, where the peculiar acceleration g is computed from 
the positions at Zi in the NAM solution. Most components agree to better than than 0.5 
percent, but there is a tail extending to —5% for one component of the initial velocity of the 
Cassiopeia dwarf spheroidal galaxy. It will be recalled that this relation assumes the particle 
orbits satisfy the equation of motion back to a(0) = 0, with coordinate displacements scaling 
as a(t). NAM produces a solution in leapfrog approximation, and the first step, from a = 
to ai = 0.1, is a long one. The small scatter in Figure [T] indicates that this first step has 
produced reasonably small — though not entirely negligible — errors in the initial condition 
in equation ([T]). That is, the approximation to the linear theory of growth of departures 
from homogeneity at large redshift is reasonably close, which ensures the wanted condition 
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that peculiar velocities at 1 + = 10 are growing in response to the peculiar gravitational 
acceleration that is causing the growth of structure. 

Numerical accuracy of the NAM solution from 1 + Zi = 10 to the present is checked by 
using initial positions and velocities derived from the solution for a numerical integration 
forward in time with 5000 steps uniformly spaced in a{t), ten times the number in the NAM 
solution. Figure [T] shows that the particle coordinate positions at l + 2;j ~ 10 are changing in 
proportion to a{t) to a good approximation, so initial conditions for the forward numerical 
integration are well approximated as 



Xi+ X2 
X3/2 = ^ , 



V3/2 = 03/2^3/2 



X2 - Xi 

a2 — ai 



(2) 



The positions xi and X2 are from the first two steps in the NAM solution, the initial value 
of the expansion parameter is 03/2 = (cti + ct2)/2, and 03/2 at 03/2 is computed from the 
Friedmann equation. The largest differences between present positions and velocities in the 
forward integration and the NAM solution are 1.5 kpc, for NGC 6822, and 1.0 km s~^, for 
LMC These are the orbits with the largest curvatures, at close passages of MW. The other 
differences of distances and redshifts generally are less than 10 pc and 0.01 km s~^. That is, 
we expect no problems with numerical errors. 




-40 -20 
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-30369 
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Fig. 2. — Differences between model and catalog LG galaxy redshifts and distances. The differences 
of distances in the right-hand panel are divided by the catalog measurement uncertainties. 
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3. Numerical Results 

Table 1 lists the differences between LG galaxy model and measured redshifts, lS.cz = 
c^mod — czcat, 8Jid distances, AD = Z^mod — -Dcat- Figure |2] shows the distributions of these 
differences. The model redshifts of 23 of the 27 galaxies are within the goal of 10 km s~^ 
differences from the catalog redshifts. The four exceptions — Cetus, Tucana, DDO 210, and 
the Sagittarius dwarf irregular — have model redshifts that are too small by ~ 30 km s^^; 
they make up the left-hand tail of the redshift error distribution in Figure |2} The model 
distances of these four misfit galaxies all are too large, by (Z^mod ~ -Dcat)/o"D ~ 6 to 9 (where 
ao is the catalog distance uncertainty, not the larger allowance in our measure of fit). 
These are the four largest distance discrepancies; they dominate the right-hand tail in the 
distance error distribution in Figure |2} The next two greatest discrepancies relative to the 
catalog errors are |-Dmod — -Dcat|/o"_D — 5 for ICIO and M33. The next largest are —4.5 for 
NGC 185, 3.8 for NGC 3109, and 3.6 for Sextans B. The rest are within three times the 
catalog uncertainty. 

The offsets D± of the particles from the catalog positions transverse to the line of sight 
are listed in Table 1. The offset of LMC is D_i = 3 kpc, which seems acceptable, and the 
rest are well within what seemed to be optimistic goals when the computation was planned. 
The last column in Table 1 is the velocities of the LG galaxies at 1 + = 10. They are 
comparable to reasonable-looking escape velocities from small galaxies, and arguably in line 
with the near close packing of the galaxies at Zi. 

The conversion from model galactocentric to heliocentric velocities depends on the cir- 
cular velocity of the local standard of rest. The minimization brings the circular velocity 
to 

Vc = 255.7 km s"\ (3) 

Table 2 compares measured and model proper motions of LMC, M33 and ICIO. The 
largest discrepancy is the motion of ICIO in declination, at three times the measurement 
uncertainty. This is arguably acceptable within the approximations of the model, and con- 
sidering the difficulty of the measurement. Table 4 lists heliocentric velocities normal to 
the line of sight in the directions of increasing right ascension and declination, and the 
corresponding angular velocities, for all the LG galaxies. 

The mass-to-light ratios m/Lx in Table 4 are based on the 2MASS K-band absolute 
magnitudes Lk in the Local Universe Catalog. These magnitudes are reliable for the larger 
galaxies, quite uncertain for the least luminous ones. The smaller galaxies were arbitrarily 
assigned nominal masses for the computation of derived from ui/Ly, = 50. Where that 
mass is too small to affect the solution the final mass nevertheless is changed by the insistence 
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Table 4: Masses and Proper Motions 
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PegDIG 


7.72 


56 


137.4 


-80.2 


0, 


,033 


-0.019 


21 


DDO210 


7.70 


54 


-21.5 


-267.0 


-0, 


,003 


-0.040 


22 


WLM 


8.98 


35 


204.0 


-210.5 


0, 


,044 


-0.046 


23 


SagDIG 


7.68 


52 


-114.6 


-285.5 


-0, 


,017 


-0.043 


24 


N3109 


9.51 


90 


13.3 


-117.4 


0, 


,002 


-0.015 


25 


Antila 


6.64 


50 


20.9 


-97.0 


0, 


,003 


-0.015 


26 


U4879 


8.24 
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-20.8 


-333.7 


-0, 


,003 
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27 


Sex B 


9.14 


51 


-29.6 


-290.1 


-0, 


,004 


-0.036 


28 


Sex A 


9.11 


49 


-35.4 


-251.8 


-0. 


,005 


-0.034 


'^Solar masses 


^km s^"*^ 


'^mas yr" 


-1 











of the computation on small but nonzero trial mass shifts in the approach to a minimum of 
X^, but the change is small. The ~ 19 final values m/L]^ ~ 50 signify masses that seem to 
be too small to matter. The possible significance of the curious values of some of the other 
masses is considered in Section HI 
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Fig. 3. — Model orbits in comoving supergalactic coordinates. Tlie numbers are keyed to the 
names in Tables 1, 3 and 4. The blue curves show the external actors. The dashed curves for the 
four misfit galaxies illustrate their apparently common origin. 

Figures |3] to |5] show orthogonal views of the model orbits plotted in comoving super- 
galactic coordinates. The center of mass of the LG galaxies plus the external actors is at rest 
and the origin is at the present position of MW. The labels near initial positions at expansion 
factor 1 + Zi = 10 are keyed to names in the first two columns in Tables 1, 3, and 4. The 
orbits of the two dominant Local Group galaxies, MW and M31, are plotted in red, and the 
four external actors that are supposed to give a phenomenological description of the effect 
of mass outside 1.5 kpc are shown in blue. The orbits of the four misfit galaxies with the 
poorest fits to catalog redshifts and distances are plotted as the dashed black curves, and 
the orbits of the other LG galaxies are shown as the solid black curves. 
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Fig. 4. — The same as Fig. [s] for the ortliogonal view toward SGL = 90°, SGB = 0. 

4. Discussion 

The focus of this dynamical study has been on fitting the model to measured positions 
and velocities of the 27 galaxies at distances less than 1.5 Mpc. We emphasize again that 
the measures of the external actors in Table 3 are chosen for a phenomenological description 
of the effect of external mass on the motions of the LG galaxies implied by the pattern 
of LG velocities and positions; the parameters of the external actors should not be taken 
as meaningful estimates of the properties of the named systems. We do take it to be an 
interesting project to attempt to interpret the large distance offsets normal to the line of 
sight for the actors named Maffei and Centaurus, and their curiously small masses, as an 
indication of the nature of the external mass distribution, but that is not attempted here. 

Since the LG galaxy masses are left nearly free to aid the fit we count as free parameters 
the ten LG galaxy masses that seem to be large enough to matter. The redshifts and distances 
of the external actors are close to their assigned values, meaning they do not seem to matter. 



-2 2 

SGY. Mpc 

Fig. 5. — The same as Fig. [3] for the orthogonal view from SGL = 0, SGB = 0. 

but we count as significant parameters the four perpendicular offsets and four masses of the 
actors. This adds up to 18 free parameters that have yielded tolerably good fits to six 
components of proper motion, 23 redshifts, and 18 distances. Whether the four redshifts off 
by —30 km s~^ and nine distances off by about 400 kpc are to be counted as modest succeses 
or serious failures is a matter of judgement. To the count of successful fits to constraints 
should be added the 28 x 3 conditions on the components of the initial peculiar velocities in 
equation ([T]), as illustrated in Figure [l| and the 28 initial velocities in the eighth column in 
Table 1, which are judged to be reasonably consistent with the close packing of protogalaxies 
at 1 + = 10. The count of constraints on initial conditions is not well defined — we can 
only offer estimates of reasonable initial velocities, and equation ([T]) is an approximation 
— but it adds to the case that the model is a useful approximation, albeit one that needs 
improvement. 

Of the 13 departures from tolerable fits to redshifts and distances, eight belong to four 
misfit galaxies, Cetus, Tucana, DDO 210, and the Sagittarius dwarf irregular. They have 
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similar discrepancies in redshifts and distances, and they have similar orbits (plotted as the 
dashed lines in Figs. [3] to |5]) that emanate from SGL ~ 200°, SGB ~ 30°. It may be 
significant that this is in the direction of the Local Void (Tully et al. 2008). The common 
features — low redshifts, large distances, and similar orbits — argue against the idea that 
measurement errors are to blame, and invites the speculation that they are manifestations 
of an inadequate external mass model that would have to be particularly serious near these 
four galaxies. If an adjustment of our phenomenological representation of the external mass 
could reduce the peculiar gravitational acceleration toward MW near the four misfits it would 
allow larger redshifts at lower present distances, in the direction wanted to improve the fit. 
A search for a fifth external mass capable of producing this effect has not yielded anything 
promising, however. An explanation of the enigmatic properties of these four misfit galaxies 
remains an interesting open issue. 

After the four misfits the next two largest numerical anomalies are the differences be- 
tween model and catalog distances relative to catalog uncertainties for M33 and ICIO. It is 
not unreasonable to expect that some distance measurements have errors well outside the 
catalog uncertainties, and the low galactic latitude of ICIO may make this distance mea- 
surement particularly problematic, though both Sanna et al (2008) and Kim et al. (2009) 
find distances to ICIO that are short of the model value in Table 1. The orbits of M33 
and ICIO are strongly constrained by measurements of their proper motions. These proper 
motions require corrections for the motions of the masers relative to the host galaxies. We 
can suggest no problem with this carefully done correction, but additional considerations of 
these valuable constraints on the orbits of M33 and ICIO and on the dynamics of the Local 
Group certainly will be welcome. 

Although the orbit of M33 fits demanding constraints from its redshift and proper 
motion as well as the initial conditions, the result is problematic. The weight of the evidence 
informing the Local Universe Catalog is that M33 is more distant than M31, but the model 
puts M33 closer. Also, the model has M33 approaching M31 for the first time since high 
redshift, while the evidence is that these two galaxies passed close enough to have disturbed 
the HI disk of M33 and drawn out stellar and HI streams possibly connecting M33 to M31 
(Rogstad et al. 1976; McConnachie et al. 2009; Putman et al. 2009; Richardson et al. 
2011). To be considered is that some galaxies that are quite isolated from visible companions 
nevertheless have apparently disturbed HI disks (Kreckel et al. 2011), and that models for 
stellar halos can produce stellar streams that originated at high redshift (Cooper et al. 2010; 
Wang, Peebles, and Nusser 2011), when the galaxies were closer together. Separating these 
phenomena from the possible effect of a more recent interaction of M31 and M33 is an 
interesting challenge for future work. 
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The model masses have problematic features, including the exceptionally large mass- 
to-light ratio of NGC 185, and, since the model would have similar circular velocities in 
MW and M31, the dissimilar model masses. Assessment of issues with masses must take 
account of three points. First, we have useful statistical measures of the relation between 
total galaxy masses and the stellar masses and velocities, but the scatter in the relation for 
individual galaxies is not at all well measured: it is to be explored by dynamical analyses 
such as this. Second, our measure of the fit of the model to the constraints uses nominal 
values of the galaxy masses that were thought to be particularly important for the dynamics. 
The relaxation to a minimum of the measure of fit can end up at a local minimum closer 
to the nominal mass than another local minimum at a more accurate mass. The effect is 
real even though the penalty of a mass adjustment is low. An example is the small model 
mass of LGS 3, only 10% of the stellar mass found by Hidalgo et al. (2011). This may be a 
result of the nominal choice of a mass small enough not to matter. Another example is the 
choice of nominal masses of MW and M31, 1.5 x lO^^m0 and 3 x lO^^m0. This approximates 
conventional wisdom, but the choice may have led to a local minimum of with similar 
model masses, 1.6 x lO^^m© for MW and 2.4 x lO^^mQ for M31. An attempt to adjust the 
solution to bring these two masses closer by adding to ^ penalty for a significant mass 
difference had little effect. This might be expected because changing the masses of these two 
galaxies requires consistent adjustments of the redshifts and distances of many LG galaxies, 
a slow operation by the present numerical method. The third point to consider is that our 
analysis has allowed masses considerable freedom to float to aid the fit to distances and 
velocities. This means that some erroneous choices of orbits may have been be made to fit 
the measured redshifts and distances by the choice of erroneous masses. Further investigation 
of the last two issues will require reconstruction of the model, which we may hope will be 
aided by future still tighter constraints on distances and proper motions that reduce the 
chance of including erroneous orbits. 

The fit to constraints on LG galaxy positions and velocities relaxes the MW circular 
velocity to Vc = 256 km s~^. This is 36 km s~^ larger than a conventional value and 26 
km s~^ larger the nominal central value used in the computation of x^? but it is consistent 
with the Reid et al. (2009a) measurement, 254 ± 16 km s~^. We cannot assign a meaningful 
uncertainty to this or the other numerical results from the dynamical model solution because 
our measure of fit is dominated by a few large discrepancies, and some (such as the four 
misfit galaxies in Figs. [3] to |5| certainly represent systematic errors in the model. That is, 
our measure has no formal significance. We can point out that the model value was 
Vc = 241 km s~^ in the solution with 10 LG galaxies (plus the 4 external actors) and Vc 
increased as more distant galaxies were added to the solution. That is, the larger value of Vc 
in the final model changed the relation between heliocentric and galactocentric redshifts in 




Fig. 6. — Orbits relative to MW, at the red circle, of the five closest neighbors. 



the direction wanted to aid the fit to the redshifts of more distant LG galaxies. This means 
that if Vc proved to be close to the standard value, 220 km s~^, the challenge to the model 
likely would not be related to the problematic orbits of the nearer galaxies M33 and ICIO. 

The computation allows the present positions of the particles to be offset from the 
observed positions of the galaxies on two considerations, that the approximate model should 
be allowed some error in this direction and that this affords an indication of whether the 
dark matter around a galaxy may be significantly displaced from the stars. The generally 
small values of the offsets D_l in Table 1 mean we have no evidence that the stars are not 
well centered on the dark matter, but tighter constraints on the model errors are needed for 
a significant exploration of this important issue. 

Figure [6] shows the evolution of positions relative to MW of the five galaxies now closest 
to MW. The coordinates are galactic and the lengths are physical. Positions at redshift z = 1 
are near the galaxy labels and present positions are at the filled black circles. None of the 
five has completed very much more than one orbit. The second and third columns in Table 5 
list the maximum distances rfmax from MW and the expansion parameters Omax (referred 
to a = 1 at 2 = 0) at dmax for the five LG galaxies in Figure [6] The next two columns 
list the minimum distances of these galaxies at a > Cmax and the expansion parameter at 
this minimum distance. The model orbits of the nearest satellites of M31 relative to its 
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Table 5: Nearer MW Neighbors 



Name 


""max 


'^max 


d ■ ^ 




LMC 


530 


0.59 


56 


1.00 


Leo I 


680 


0.43 


200 


0.93 


Leo T 


580 


0.71 


420 


1.00 


Phoenix 


760 


0.60 


460 


1.00 


N 6822 


610 


0.43 


34 


0.79 


""kpc 



position have similar appearance to Figure [6} though none pass as close to M31 as do LMC 
and NGC 6822 to MW. This should be taken with caution, however, because the present 
distances from MW to its nearest neighbors are reasonably well known, while modest relative 
distance errors could introduce considerable errors in the present distances from M31 of its 
nearer neighbors. 

In the model, the galaxies Leo T and Phoenix are approaching MW for the first time and 
LMC is now passing MW. The galaxy Leo I is moving away after closest approach of 200 kpc 
at redshift z = 0.08. Leo I has been compared to closer dwarf spheroidal satellites of MW 
that have completed many orbits (and for this reason are not in the compuation) and in the 
process may have been expected to have lost much of their gas. Leo I does not have detected 
HI (Young 2000), despite the single relatively distant passage at 200 kpc distance. If the 
model orbit is right it means either that the ram pressure of plasma clouds at this distance 
was capable of removing the HI, or that Leo I was capable of driving away the remnant HI by 
itself. The predicted much closer passage of NGC 6822, 34 kpc at redshift z = 0.27, left HI 
around this galaxy, and might be expected to have left a tidal tail similar to the Magellanic 
Stream, though perhaps significantly dissipated by the passage of time. It may be significant 
that the redshift in the HI around NGC 6822 varies almost linearly with position along the 
long axis of the HI distribution, not the expected behavior of a rotationally-supported disk 
in a standard dark matter halo (Weldrake, de Blok, & Walter 2003). To be investigated is 
whether the observed distribution and motion of the HI and stars in this galaxy might be 
consistent with remnant effects of a close MW passage. 

Table 4 lists proper motions of the LG galaxies. Van der Marel & Guhathakurta (2008) 
deduce from the motions of its satellites that the transverse components of the heliocentric 
velocity of M31 are f „ = 78 ± 41 km s~^ and vs = —38 ± 34 km s~^. The former differs 
from the dynamical model by four times the stated uncertainty and the latter agrees with 
the model within the uncertainty. Since the model has its own uncertainty, though we have 
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no meaningful estimate of it, the two measures of the hehocentric velocity of M31 seem to 
be tolerably consistent. 

Measurements of proper motions of the nearer galaxies in the Local Universe Catalog 
would be particularly important additions to the constraints on local dynamics because 
their distances from MW are fairly well known, and their model orbits reach far enough 
from MW to be to be sensitive to the evolution of the local mass distribution on scales of a 
few megaparsecs. The dwarf elliptical Leo I may be an interesting candidate for an optical 
measurement of the mean angular velocity of the stars. In the model it has the second 
largest heliocentric angular velocity, 0.33 mas yr~^, about three times the uncertainty in the 
measurement of the angular velocity of LMC Half of the angular velocities of LG galaxies 
in Table 4 are greater than 0.05 mas yr~^, well above the detected motions of masers in 
ICIO and M33 (Brunthaler et al. 2005, 2007). Depending on what searches for masers in 
LG galaxies reveal, measurements of angular motions by present VLBI technology has the 
potential to open a new frontier in constraints on local dynamics. We may hope this frontier 
will be opened far more widely by new advances in the optical (the GAIA Science Mission) 
and VLBI (Reid et al. 2009b). 
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Fig. 7. — Positions of the LG galaxies at 1 + = 10 in physical supergalactic coordinates. The 
ten galaxies in the concentration at SGL~ 70°, SGB~ are plotted as filled squares. 



A group of ten LG galaxies emanate from high redshift along similar paths. This is 
seen in Figures [3] to [5} and it is illustrated in another way in Figure [7], which shows physical 
positions at the starting time of the computation, at expansion factor 1 + = 10. The 
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external actors are not shown in this figure, the initial positions of MW and M31 are plotted 
in red, and the LG galaxies not in the group of ten are shown as filled circles. The ten with 
similar initial positions, shown as filled squares, are in order of increasing present distance 
from MW the galaxies NGC 6822, NGC 185, NGC 147, Leo A, Andromeda XIV, Cassiopeia 
dwarf spheroidal, Pegasis dwarf, UGC 4879, and Sextans A and B. One might include a more 
distant eleventh member of the group, Andromeda XII. It is plotted as the filled circle at 
largest SGX and SGY in Figure [7| and has the label 17 in Figs. [3]to|5| The initial positions 
of the inner ten are in a band of physical length ~ 400 kpc that is close to perpendicular 
to the supergalactic plane, and centered at SGL ~ 70°, SGB~ 0. The narrow dimension 
is less than 75 kpc. If this distinct early concentration had been a little tighter the group 
could have merged into a protogalaxy. It might prove interesting to check whether these ten 
parts of a possible failed protogalaxy have more in common among themselves than among 
the other low mass LG members, perhaps in the HI mass fraction or the distributions of 
metallicity or stellar evolution ages, refiecting a near common origin. 



5. Concluding Remarks 

This model for the Local Group certainly has problems. Perhaps most serious is the 
orbit of M33, which fits the measured distance, redshift and proper motion, within tolerable 
errors, along with the initial condition, but not the indications that M33 passed close to 
M31 at modest redshift. This situation certainly requires further examination. The common 
anomalies in the redshifts and distances of the four misfits whose near common orbits are 
plotted as dashed lines in Figures [3] to [5] seem to be pointing to an error in the simplified 
phenomenological description of the effect of the mass distribution outside the Local Group. 
Investigating that will require a more detailed dynamical analysis that takes explicit account 
of where the larger external masses are observed to be now and where they are computed to 
have been. We hope to report on an exploration of this analysis in due course. The values 
of some of the masses derived from the dynamics are curious, and may be in part the result 
of compensation for systematic errors in the model. 

The very real problems should not obscure the evidence that this dynamical model has 
proved to be a useful approximation to reality, from the success in fitting the considerable 
number of observational constraints on positions and velocities in Tables 1 and 2, and perhaps 
also the Milky Way circular velocity Vc in equation (|3|, with orbits that have the initial 
conditions required by a reasonable approximation to the cosmological growth of structure. It 
is notable also that the model has led to the considerable list of issues for further consideration 
presented in Section |4j 
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In the past decade advances in detections of nearby galaxies, and in the measurements 
of galaxy distances and velocities, have greatly enriched the study of dynamics in the now 
densely sampled region within 1.5 Mpc distance. We have the prospect of another advance 
of the same order from work in progress, or that seems to be feasible, on measurements 
of galaxy distances and proper motions. That will present us with a really serious but 
fascinating challenge for the development of numerical methods of analysis of a tangled web 
of constraints on what happened in the local universe. 



We are grateful to Mark Reid and Jeremy Darling for guidance to prospects for mea- 
surements of proper motions of Local Group galaxies. 
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A. Yet another Version of the Numerical Action Method 

This version produces a solution to the equation of motion in leapfrog approximation 
by the joint relaxation of all coordinatess toward a stationary point of the action in the 
direction indicated by the first and second derivatives of the action. This joint relaxation of 
all coordinates describing all orbits considerably speeds the computation when applied to a 
modest number of particles. 

In the cosmologically fiat universe of this analysis the expansion parameter satisfies 

"''' + (i-mi ^ = -^+(i-fi)//j, (Ai) 



with present value Go = 1. 

Particle orbits are represented by comoving positions Xi^k,n for the particle label 1 < 
i < Up, cartesian component 1 < k < 3, and time step n < rix+i- The earliest computed 
positions are at ai = Oj, the given present positions are at a^+i = 1, and the coordinates to 
be adjusted are at 1 < n < n^^. The leapfrog numerical integration commences at ai/2 = 
at the nominal singular start of expansion. 

The equation of motion in leapfrog approximation is 

^ ^ _ Qn+l/2Qn+l/2 _ 



+ 



tn+1/2 — tn-1/2 



(A2) 



an-l/2«n-l/2 _ 

' \''^i,k,n •'^i,k,n—lj- 

^n— 1 



The expansion time tn+1/2 is computed from the analytic expression for the time from qq — 
to a„+i/2 = (on + an+i)/2. The physical acceleration of particle i produced by the gravita- 
tional attraction of the other particles is Qi,k,n/o^i- The physical acceleration corresponding 
to the counter term in the square brackets may be written as 

Gpnri,k,n, (A3) 



2al 3 

where r = ax is the physical coordinate and p„ is the cosmic mean mass density at a„. The 
counter term thus causes the peculiar acceleration to vanish when the physical acceleration 
matches that of a homogeneous universe. When the particles are represented as point masses, 
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In the limiting isothermal sphere model for the mass distribution in MW, i = 1, when there 



is a nearby smaller galaxy, j, the terms for this pair in the sums in equation (A4) are replaced 
with 



8Qjh,n = Grrii- 



(A5) 



at separation 

\xi,n - Xj^nl < Xo = Gnii/ {vlan) ■ (A6) 

This means the physical acceleration of particle j caused by MW is f^/r, and the acceleration 
of MW caused by particle j is irtj/mi times v1/r, which conserves momentum. 



Equation (A2) is simplified by setting 



^n-l/2'^«-l/2 _ + dtn _ in+1/2 " ^ri-1/2 

n— 1' ■ 



where 



Ff = = Fo+ 



follows from the conditions ai/2 = and a oc 

^2/3 

at a — i- 0. This brings eq. <\A2h to 



1 

_( 

2^ 



(A7) 
(A8) 

(A9) 



The coordinates are relaxed toward a solution at Si^k,n = by the coordinate shift Sxi^k,n 
that satisfies 

Si,k,n + Si^k,n;j,k',n'SXj^k',n' = 0. (AlO) 

j,k' ,n' 

In this expression the nonzero derivatives of the Si^k,n with respect to the coordinates are 



^i,k,n;i,k,n+l 
Si^k,n;j,k' ,n 

Si k,n:i,k' ,n 



i,k,n:i,k,n—l 



Gi,k,n;j,k', for j ^i, 



(F+ + F-)5k,k' + 



dtr 



Gi,k,n;i,k' + ^^H'^Sk,k' 



The derivatives of the acceleration for i ^ j are 



Gi,k,n;j,k' — GlTlj 



for the inverse square law, and, for the isothermal sphere model. 



Gi,k,n;j,k' — Guij 



(All) 



(A12) 



(A13) 



For i = j the derivatives are 



Gi,k,n;i,k' ^ ^ Gi,k,n;j,. 



Uk' 



Equation (AlO) with the nonzero terms ehminated is 



Setting n + 1 — )• n and rearranging gives 

Si^k,n~l ~l~ X^^ifc' ^i,k,n—l;j,k',n—l^-^j,k',n—l ~l~ Si^}^^Yi—l;i,k,n—2^-^'i,k,n—2 



^•^i,k,n 



which iterates to 



At n = 1, 



Si^k,n—l\i,k,n 

SXi^k,n = Ai^k,n + Bi^k,n;j',k"SXj'^k",l- 
j',k" 



AihA — 0, BikA-j'h" — ^Lj'S. 



j"Jk,k" 



At 77. = 2, equation (A16) is 



Si,k,l + Si^k,l]j',k",l^Xj'^k",l I Si^k,l]i,k,2 

j',k" 



SO 



Ai,k,2 — —Si^k,l/ Si^k,l;i,k,2, -Si,A:,2;j,fc" — ~ Si^k,l;j,k" ,1 / Si^k,l■,i,k;. 



At Tl > 3 the form (|A17|) in equation (|A16|) gives 



(A14) 



Si^k.n ~l~ ^i,k,n;i,k,n+l^Xi^if^fi-\-l ~\~ ^ ^ Si j^ n-j k' Yi^Xj j^i fi -\- Si k,n—l^X.i j^ Yi—l 0. (Al5) 



(A16) 



(A17) 



(A18) 



(A19) 



(A20) 



Si±,n-l + 'S'j^fc^„_i;j_fc/^„,_i ^Aj,fc',n-1 + Bj^k/^n~~l;j',k"^Xj'^k",l 

j,k' j',k" 



(A21) 



Thus at 3 < n < the coefficients are 

Si^k,n—l ~l~ X!]ifc' '5'j,fc,n— l;j,A;',n— l^j,A;',n,— 1 ~l~ Si^k,n—l;i,k,n—2-^i,k,n—2 



"-"ij/c,?!— l;j,fc,n 

Xlifc' Si^k,n-l;j,k' ,n-lBj^k' ,n~l;j' ,k" + Si^k,n-l;i,k,n-2Bi^k,n-2;j' ,k" 



i,k,n;j' ,k" 



Si k,n- 



l:i,k,n 



. (A22) 



This gives the coefficients v4j_fc „ and Bi^k^n;j' ,k" by iteration starting from equations (A18) 
and (|A20|). 



Equation (A15) at n = n^, with Sxi^k,n^+i = because present positions are fixed, is 



or 



j,k' 

= Si^k,nx "I" ^ ^ Si^k,nx;j,k' ,nx ^j,k' .Ux ~^ ^ ^ Bj^ki^nx;,j',k"^Xji^k",l 



^i,k,nx~l + Bi^k,nx-l;j',k"SXj'^k",l 
j',k" 



— Tj^fe + Ti^k;j',k"SXj'^k". 



1) 



-^j,fe Si^k^nx ~l~ ^ ^ Si^k^nx;j,k' ,nx-'^j,k' ,nx ~^ Si^ji;^nx',i,k,nx — l-'^i,k,nx — ly 



j,k' 



Ti^k;j',k" ^ ^ Si,k,nx;j,k' ,nx-Bj^k' ,nx;j' ,k" ~l~ Si ]f Yix\i,k,nx — l-Bi^k,nx — l;j',k"- 
hk' 



(A23) 



(A24) 



Inversion of the matrix Ti^k;j',k" fixes the ^Xj^fc^i, and equations (A17) and (A22) give the rest 
of the Sxi^k,n- 



